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We consider the finite temperature Casimir eff'ect between two concentric spheres due to the 
vacuum fluctuations of the electromagnetic field in the {D + l)-dimensional Minkowski spacetime. 
Different combinations of perfectly conducting and infinitely permeable boundary conditions are 
imposed on the spheres. The asymptotic expansions of the Casimir free energies when the dimen- 
sionless parameter e, the ratio of the distance between the spheres to the radius of the smaller 
sphere, is small are derived in both the high temperature region and the low temperature region. It 
is shown that the leading terms agree with those obtained using the proximity force approximation, 
which are of order Te^^^ in the high temperature region and of order e^^ in the low temperature 
' region. Some universal structures are observed in the next two correction terms. The leading terms 

of the thermal corrections in the low temperature region are also derived. They are found to be 
(~| , finite when e — >■ O'*', and are of order T^^^. 
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^ \ I. INTRODUCTION 

Ph. 



Casimir effect is one of the most interesting macroscopic phenomena in the quantum theory of fields. It has been 
under active studies under various context The success in its experimental verification |2t21 has intensified the 
^2 , interest in this effect. In recent years, the thermal correction to the Casimir effect has attracted increasing interest 
Q^' both theoretically and experimentally P, 0-[i3 • 

The interest in the Casimir effect of spherical objects can be dated back to the work of Boyer [ll| , where he computed 
■ the zero temperature Casimir force acting on a three-dimensional perfectly conducting spherical shell and found that 
^1 it is repulsive. This result has later been confirmed in a number of other works |12h19I |. Since then, the Casimir 
• effect in spherical configuration has attracted considerable interest. The eases of scalar fields, spinor fields and vector 
' fields have been considered in various works pol - l28j . In fact, in (2TI426I . [28j. the authors considered spherical shells 
. in general {D + l)-dimensional Minkowski spacetime rather than restricted to 4-dimensions. One of the motivations 
for this is that physics in higher-dimensional spacetimes have become a trend since it was found that the existence 
, of extra dimensions might be used to resolve some fundamental problems in physics such as the hierarchy problem. 
' Another reason is that the dimension of spacetime can be used as a perturbation parameter in quantum field theory 
calculations [2l|,[2l[3^. 

For the last ten years, there has been an intense interest in studying the Casimir interaction between two objects. 
Several methods have been developed to compute the Casimir interaction beyond the proximity force approximation, 
such as the functional determinant or the multiple scattering method (3ll-l36j and the worldline approach [37i - [39| . 
The corrections to the proximity force approximations have been computed for several geometric configurations such 
', as the sphere-plane [3l|, [H, liol - lisj , cylinder-plane [3l|, HI, |43 - I46j , cylinder-cylinder [46l - l48j , sphere-sphere [i^ [s^] , 
etc. Recently, there has been an interest in considering the finite temperature correction to the Casimir interactions 

PH. 

As a matter of fact, the Casimir interaction between two concentric spheres in D = 3 dimensions has been con- 
sidered in (56l - [5^ . For scalar fields or spinor fields in general D-dimensions, the zero temperature Casimir effect on 
two concentric spherical shells has been studied in [sqHb^ . In [65j . we have derived the zero temperature Casimir 
interaction between two concentric spheres due to the fluctuations of electromagnetic field in the D-dimensional space. 
In this article, we consider the finite temperature effect. Moreover, we would derive the asymptotic behaviors of the 
Casimir free energy when the separation between the spheres is small. The results are compared to the proximity 
force approximations. 

In this article, we use units where h = c = = 1. 
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II. CASIMIR FREE ENERGY OF CONCENTRIC SPHERES 



Consider two concentric spheres with radii ai and 02 (oi < 02) in (D + l)-dinicnsional Minkowski spacetime, 
with either the perfectly conducting or the infinitely permeable boundary conditions. The electromagnetic field 
Ffj,u — d^A^, — d^Afj, satisfies the field equation: 



-^5^ {V\9\9''^9''\dxA^ - dM) = 0. 



As usual, the Coulomb gauge 



= 



is imposed to remove the gauge degree of freedom. The eigcnmodes of the field arc divided into TE modes and TM 
modes. In terms of _D-dimensional spherical coordinates (r, 6),9 ^ {9i, . . . , (?d-2, 0), the TM modes have the form 



dr 



The TE modes can be divided into {D — 2) sets, each of them has the form 

Ar =0, 

D — 2 

Here v ~ I -\ — , m = (to2, • • ■ , niD-i), Jv{z) and N^{z) are Bessel functions of first kind and second kind. The 

sets of m have been discussed in detail in [6^1 • For TM modes, each fixed I has 

{2l + D-2){l + D-i)\ 
^ {D^l 

allowable m; whereas for TE modes, each fixed I has 

, l{l + D~2){2l + D-2){l + D-4)l 

^'^^^ = (z?-3)!(; + i)! 

allowable m. The perfectly conducting boundary condition is equivalent to 

'dA„, dAg 



-0, 



l<i<i <D-l. 



boundary 

Therefore, for TE modes, the perfectly conducting boundary condition on r = implies 

CiJ^{a,uj) + C2N^{a,uj) = 0, 

whereas for TM modes, we have 



0. 



The infinitely permeable boundary condition is equivalent to 



fdAe, 


dAr\ 


\ dr 





0, 



l<i< D~l. 



boundary 

Therefore, for TE modes the infinitely permeable boundary condition on r = a; implies 



dr 



(r^ [CiMruj) + C2N,{ruj)]) 



0, 
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whereas for TM modes, we have 



CiJ^{a,uj) + C2N^{a,Lo) = 0. 

The interacting Casimir free energy E'cas of the concentric spheres can be written as a sum of the contribution from 
the TE and the TM modes: 



Cas — ^'TE + ^'TM 



(1) 



In the foUowing, we wih use E to represent either i^xE or -Etm- Using zeta regularization method, the Casimir free 
energy is given by (see e.g. eq. (5.17)): 

i^Ca. = -| (Ct(0) + log[/i']CT(0)) , (2) 

where ^ is a regularization parameter with the dimension of mass, Ct('5) is the zcta function 

oo 

UJ P— — 00 

uj are the TE or TM eigenfrequencies, and = 2ttpT are the Matsubara frequencies. As in the zero temperature case 
iGal. one can show that 



(3) 



(4) 



Q • poo 1 

= y: md) E' / ie - ep)-'^ji{Od^, 

1=1 p=0 ''^p ^ 

where di{D) = hi{D) for TE modes and di{D) = hi{D) for TM modes, and 

/KC)=ln(l-A/KC)), 

/i/(z) and Ku{z) arc modified Bessel functions of first kind and second kind. The values of Ui and j3i depend on the 
type of modes and the boundary conditions imposed on the sphere r = ai. They are hsted in Table HI 



TABLE I: The values of ai and jSi under different boundary conditions. 



Type of mode 


Sphere r = ai 


ai 


ft 


TE 


perfectly conducting 


1 





TM 


perfectly conducting 


D-2 
2 


1 


TE 


infinitely permeable 


A- D 
2 


1 


TM 


infinitely permeable 


1 






As in [65|, one can show that 



1 = 1 p=o ''ip 

is an analytic function of s. Therefore, Ct(0) = and 



oo OO „oQ J OO OO 

c^(o) = 2Erf;(^)E' / ^'(^)'^^ = -2E'^'(^) E/'(^p)- 

1=1 p=o ''ip ^ 1=1 p=o 
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Since 

as z — >■ 0, we find that 

/i(0)=ln 

Therefore, from ([2]), we find tliat tlic TE or TM contribution to tlie Casimir free energy is given by 





+ I3iv){a2 






2v\ 




- I3iv){a2 




.«2. 





2 











[ai - /3ii^)(q;2 









(5) 



=1 p=i 



Using Poisson summation formula, the Casimir free energy ([5]) can also be written as 

1=1 p=l -^0 ^ 



(6) 



where 



1 _ _ 

is the zero temperature Casimir energy (vacuum energy). The expression ([5]) is suitable for the study of the high 
temperature limits of the Casimir free energies, but for the low temperature limits, the expression (|6]) would be 
preferred. 

From the expression ([5]) for the Casimir free energy, one can use the argument in [65| to show that the force acting 
on the spheres is always attractive when the two spheres have the same boundary conditions (homogeneous boundary 
conditions); and is repulsive when one of the spheres is perfectly conducting and the other is infinitely permeable 
(mixed boundary conditions). 

III. PROXIMITY FORCE APPROXIMATION 

In this section, we discuss briefly the proximity force approximation of the Casimir free energy when the separation 
between the spheres is small compared to both radii of the spheres. Define a dimensionless parameter 

a2 — ai d 
ai ai ' 

where c? = 02 — ai is the distance between the two spheres. In the following, we are going to study the asymptotic 
behaviors of the Casimir free energy when e ^ 1. We consider the following two regions: 

1. Low temperature: dT <ti a2T <C 1; 

2. High temperature: a2T > dT > 1. 

For a pair of infinite parallel plates in {D + l)-dimensional spacetime, the Casimir free energy density, as a function 
of the separation between the plates d, is given by [g^, 113] : 



sUd) ^{D-l)\ -^^^^UD + 1)^ + ■■■). (8) 



r(^) 1 
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if both the two plates are perfectly conducting or infinitely permeable. We refer to these as the homogeneous boundary 
conditions. ([5]) is the low temperature asymptotics, the leading term being the zero temperature term; and (jH]) is 
the high temperature asymptotics, the leading term being called the classical term. The factor {D — 1) is due to the 
{D — 1) polarizations of photons in {D + l)-dimensional spacetime, (D — 2) of them come from the TE modes and 
one of them comes from the TM modes. 

In case of mixed boundary conditions, i.e., one plate is perfectly conducting and one plate is infinitely permeable, 
the corresponding Casimir free energy density is 

4as W = - 1) - 2-^Kr{D + 1)^ + . . . j , 

or 



In case of two concentric spheres, proximity force approximation of the Casimir free energy is particularly simple. 
It is given by the product of the surface area of the sphere (either one) with the Casimir free energy density between 
two parallel plates. Since the surface area of a (-D — l)-dimensional sphere of radius ai is 



proximity force theorem implies that in the low temperature region, the proximity force approximation to the Casimir 
free energy between the spheres is 



for homogeneous boundary conditions, and is 

Ellt - El--, = ^T^l^^d - ^-""KniD + 1)^^ (11) 

for mixed boundary conditions. Notice that these leading terms are of order e^^, D being the space dimension. 
In the high temperature region, the proximity force approximation gives 

Ell^ - Ei:^-' = -^j^T^Cn{D)T (12) 
for homogeneous boundary conditions, and 

Ell^ - E'cl^-' = ^^^TT^^H 2^-^)U{D)T (13) 

for mixed boundary conditions. These are of order Te^~^ . 

In the following, we use the exact formulas ([5]) and ([6]) to find the asymptotic expansions of the Casimir free 
energy and compare to the proximity force approximations (fTOl) . (fTT|) . ([T2|) . (fT3|) . We start with the high temperature 
asymptotic expansion because it is less technical. 
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IV. SMALL SEPARATION ASYMPTOTIC EXPANSIONS OF THE CASIMIR FREE ENERGIES IN 

THE HIGH TEMPERATURE REGION 

From ([5|), it is easy to see that in the high temperature region where ^ dT 3> 1, the Casimir free energy is 
dominated by the term (classical term): 



T 



1=1 



("1 


+ I5iv){a2 




ai' 


2iy\ 


("1 


- I3iv){a2 




.«2. 





(14) 



In the following, we derive the asymptotic expansion of this term when e <C 1 . We consider the case of homogeneous 
boundary conditions and the case of mixed boundary conditions separately. 

A. Homogeneous boundary conditions 

In the case both spheres have the same boundary conditions, we find that the high temperature limit (the classical 
term) of the Casimir free energy (|14p is the same and is given by 



Ol 

a2 



T 



-Y^di{D)\n{l 



-2au\ 



1=1 



where 



iog^ = iog(i + -)=E(-ir - 



a2 



ai 



i 



(15) 



Since hi{D) and bi{D) can be expanded as 



D-2 



hi{D) = ^ XD-jv\ 



]=0 



D-2 



D-2 



we can write di[D) as di{D) = zuu-ji/^ . Then 



3=0 

rj-^ D — 2 OO OO 

rj-i D—1 pc-{-ioo 

i=i 



r(z)(2a)-Xfl(2 + i)Ch + 1; f ) dz 



21 



We have used the inverse Mellin transform formula 



-1 /.c+ioo 

e-" = : / r(z)u-^dz, 

27rz J ^_Ar^ 



(16) 



(17) 



and the fact that the Hurwitz zeta function C,h (s; c) = ^~^(n + c) * has a single simple pole at s = 1 with residue 



n=0 



1. In fact, the use of the inverse Mellin transform p7)) and the residue theorem allows us to find the full asymptotic 
scries in a from the second line of (jl6p . The last line (|16p gives the asymptotic scries up to the term in In a. Using 
(|15p. we can rewrite this asymptotic expansion in terms of e. Since 



xd-d-2 



(£)-3)!' 
2 

yD-D~2 —71=; 7;^, VD-D-Z — 0, VD-D-i 



{D-2] 



D^-6D + 32 
12(i:>-4)! ' 
1 

12(1? -5)!' 



(18) 
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we find that when D > A, the first three leading terms of the TE and TM contributions to the classical term are given 
respectively by 

Etb-E^e ^l + e^ + e s ^^^^^^ 

pel pPFA.eiL ^^^-1 o(3£-8)CD-^ 

where E^^''^^ and i?^^''^' are respectively the proximity force approximations for the TE and TM contributions 
which are {D — 2)/{D — 1) and 1/{D — 1) times the total proximity force approximation ((T^ . The first three leading 
terms of the Casimir free energy is then given by 

cl . c.PFA,df, , D-1 , 2(3D-8)(D-1) ^D' -bD + 2SU(P-2) 



When D = 3, a more precise computation shows that the first three leading terms of the Casimir free energy is 
ES^ - bS^- { 1 + . + ^/ ,„.} _|,a(3, { 1 + . + 1„ . 

This has equal contributions from the TE and the TM modes. From (fTg]) . we see that proximity force approximation 
underestimates the Casimir free energy, and the underestimation is worse when the space dimension becomes larger. 



B. Mixed boundary conditions 

In the case r = oi is perfectly conducting and r = 02 is infinitely permeable, we find that the classical terms of the 
TE and TM contributions to the Casimir free energy are given respectively by 



4-r> 



1=1 



ai 


2i/\ 


.«2. 





1=1 



In the case r = ai is infinitely permeable and r = 02 is perfectly conducting, we have 



4-D 



1=1 



Consider series of the form 



/ = ^di(Z?)ln 1 



1=1 



4-D 





2iy^ 






fli 









D-2 



D-2 



D-2 



ai 


2i/^ 






ai 




.«2. 





which can be rewritten as the sum of two terms: 

D-2 00 , ^ 

1 = Y^ ^D:j E (mi + e"'"") + In f 1 + 
i=o /=i ^ ^ 

D — l 00 00 

= E^^u-iE^'"'E 



2i9 



1 



^) -2ak 



V --d e^"" + 1 

D-l a 



1=1 k=l 



J = l 



1=1 



h+h 



(20) 



For the /i term, we find as before 



D-l 



VUD;j-l 



r(j)(l-2-^)Cii(j + l)- + 0(lna), 
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where we have used 



fc=i 



For the I2 term, since 



00 1-1 1 poo D-2\i^^ 1 

V — ^ ^ / + — ) ^ dx 

J2 r + " M {2a^yX,.k-r-uk + 0(ln a), 



where 



we find that 



^ ' r=0 



j = l fe=l ^ r=0 \ / 

From these, we can derive the leading terms of the Casimir free energy. 

D - 4 

If the sphere r — ai is perfectly conducting and the sphere r — a2 is infinitely permeable, then i? — — - — and 
D-2 

1!} = — - — respectively for TE and TM modes. When D > 4, wc find that the first three leading terms of the TE 
and TM contributions to the classical term are given respectively by 



PFA,cil, , D-1 2{D-A)2" -%Cb{D-2) ^{ZD - %){D ~ I) 



D-2 2^-2 Cr{D) 24 

^50"^ ~30D + 16 2^ ^SCniD -2) ^ 2{D - 4^ 2^ - 32 Cr{D- 4) 
+ ^ 6{D-2) 2^-2 Cr{D) +^ {D-2)(D-3) 2D -2 Cr{D) 

Etm-E^^ Y + e^ + 2s^^^-^+e 

')D _ Q /-^f ri _ o\ of n _ o\ oD 



(21) 



^ 50-82^ ^8CRiD-2) , 2(Z3-2) 2^ - 32 Cfl(i^ - 4) 
6 2^-2 Cr{D) iD-3) 2D -2 ^r{D) j' 

where E'^^'^''^' and -E'^^''^' are respectively the proximity force approximations for the TE and TM contributions 
which are {D — 2)/{D — 1) and 1/{D — 1) times the total proximity force approximation ([T3|) . The first three leading 
terms of the total Casimir free energy is thus 

pel pPFA.ei/,^ ^-1 2p - 3) 2^ - 8 Ci^(J - 2) , ,{3D-8){D-1} 



^•^^ Cas \ 2 D-l 2^-2 Cr{D) 24 

(22) 



_2 



5L»2 _ 25D + 8 2^ - 8 CflXi? - 2) 2 2(£'^ -7D + 14) 2^ - 32 Cr{D - 4) 



6{D~1) 2D -2 Cr{D) {D-1){D^3) 2^-2 Cr{D) 

When D = 5, the term (2^ - 32)(r{D - 4) in (gT]) and ^ is understood as 

lim (2^ - 32)Cfl(i:' - 4) = 32 In 2. 

D->5 
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When D = 3; a more detail computation gives 



_2 I _ _37^ 

3^'"" 3C^ " T6i^^"^"M ^ ' ^ 3^'"" 3Cfl(3) 



^^TE -E^l^''' J 1 + e - In 2 - in e ^ = Jii^Ci? (3) <| 1 + e - |e In 2 - lne\, (23) 



I Q 2 1 I 8 2 1 



3r 



3Cfl(3) j 8£2^"^^|^ 3Cfl(3) 

4 - L» 

If i/ie sphere r — ai is infinitely permeable and the sphere r ~ a2 is perfectly conducting, then d = — - — and 
2- D 

i!} ~ — - — respectively for the TE and TM modes. When D > 4, the first three leading terms of the TE and TM 
contributions to the classical term are given respectively by 

pel pPFA,ci/, ^^^-1 J{D-4)2^-8Cr{D-2) , ^, {3D-8){D-1) 
^™-^^TE Y + '^T^' D-2 2D-2 Cr{D) +^ 24 

^2 VIJ^ - 42£) + 80 2^ - 8 Ci?(-D - 2) _^ ^2 2{D ~ A)"^ 2^ - S2 Cr{D - 4)\ 



6{D~2) 2D -2 Cb{D) {D-2){D-2,) 2D -2 Cr(D 

Etm ^ 1 + £ — 2e 



ci ^PFA,cil, , D-1 2D-8UD-2) , ,{3D-S)iD-l) 



2D -2 CbXD) 24 

2 ?£> - 16 2^ - 8 UiD - 2) ^2(0 - 2) 2-° - 32 Cr{D - 4) 1 
6 2^-2 CbW) +^ {D-3) 2D -2 Cr{D) j' 

The first three leading terms of the total Casimir free energy is thus 

pel pPFA,ci/,^ D-l_ 2{D-3)2D-8Cn{D-2) , i3D-8)iD-l) 

^Cas-^^Cas j^+S 3 " - 1 2^ - 2 Ck(^) 24 

2 TL'^ _ S5D + 64 2^ - 8 Ci?(£' - 2) 2 2(-Ci^ - TZ? + 14) 2^ - 32 Ci?(-D - 4) 
6(Z?-1) 2^-2 Cfl(^) (I?-l)(i?-3) 2^-2 CflP) 



When D = 3, by duality, the TE contribution is given by (P^ . and the TM contribution is given by 

Observe that the corrections to the proximity force approximation in the case of mixed boundary conditions is more 
complicated than the case of homogeneous boundary conditions. We also find that the proximity force approximations 
underestimate the Casimir free energies. 

Compare the first correction to the proximity force approximation for the two scenarios of mixed boundary con- 
ditions, we find that in the high temperature region, the force is stronger when the sphere with smaller radius is 
perfectly conducting and the sphere with larger radius is infinitely permeable. 



V. SMALL SEPARATION ASYMPTOTIC EXPANSIONS OF THE CASIMIR FREE ENERGIES IN 

THE LOW TEMPERATURE REGION 

In the low temperature region, the Casimir free energy is dominated by the zero temperature term. Making a 
change of variables ^ h- !■ uj/ai in ([7]), we find that the zero temperature Casimir free energy can be written as 



2^ 00 „oo 

Eo^- ydi{D) \n{l - A,{iu))dio 

27rai ^ Jo 

^ s=l 1=1 
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where 



{aih{u) + piLoIliuj)) ia2K4uj{l + e)) + + + e))) 



iaiK,{io) + piLoKliio)) (a2/.(w(l + e)) + + e)/^(w(l + e))) 

From Debyc asymptotic expansions of Besscl functions (68l . l69j , we have 



aiKu{vLo) + l3iVioKl{vijj) 



{ I ( 

— exp 

\ k=l 



,2k-l 



, if ai = 1, Pi = 

- exp I 2i'7](u}) + 2 2^ - 



,2fc-l 



if /3. = 1 



fc=i 



where 



7(2) =Vl + z2+log 



1 + 



t(z) = 



vr 



V^.mfi + V-W 



k=l \ k=l / k=l 

Uk{t) and Vfe(t) are polynomials in t defined recursively by 

t2(l-i2) 



fc=l 



uo{t) = 1, Ufc(t) = 



In the following, we discuss the asymptotic expansions of the zero temperature Casimir free energy for the case of 
homogeneous boundary conditions and the case of mixed boundary conditions separately. 



A. Homogeneous boundary conditions 

In this case, ai = a2,Pi = (32- We find that 



exp I -2i^ (7;([1 + e]u;) - 77(0;)) -^2^^ ' J^SF^T I 



(25) 



The polynomials Pk{t) are equal to Dk{t) or Mk^a{t) depending on the type of modes and the boundary conditions, 
as shown in Table HIl 



TABLE II: The polynomial Pk{t) under different boundary conditions. 



Type of mode 


Boundary conditions on both spheres 




TE 


perfectly conducting 


Dk(t) 


TM 


perfectly conducting 


M^O_2{t) 


TE 


infinitely permeable 




TM 


infinitely permeable 





In the following, we only find the first three terms in the asymptotic expansion when e — > 0+. For this, it is sufficient 
to take the /c = 1 term in (|25p . The polynomials Di{t) and Mi^a{t) are given explicitly by 



Dl(t) = , 

^' 8 24' 
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Therefore, we can write Pi{t) as Pi{t) = Xt + -ft^. Using p7|) . we have 



En 



oo D~2 
oo D~2 



J + 1 



exp —2s 



I/(7?((l+eV)-77H) + 



PlWc.(l +£)))- Pi (tH) 



1 CX3iV — ^ pOO -I nC-\-i00 



X 1 + 



(z + l)z.-^ (77((l+eH-77H)- 
1 Pi(t(L^(l+e)))-Pi(<M)~ 



v'^ ri{{l + e)u}) - ri{u}) 



dzduj. 



Since 



(?7((1 + £M -r/H)- 



e uj 



1 



1 Pi(t(c^(l + e)))-Pi(tH) 



we find that 



2 r;'H 
, z Pi'(tM) 



6 7/(w) 8 



D-2 



Eo ^ - - — ^D-j-z— 

27rai ^ 27ri 



where 

/•oo 

^(z) = / {co^iico)) 
Jo 

nOO 

B[z) = / {u^^iiu^)) 
Jo 

It is straightforward to find that 

A{z) 



^_ euj r]"{uj) 

^ 2 7?'(C.) 

P{(t(c.))t'(.. 



za;2 ?7"'(ii^) z{z + 1)0;^ ?7"(w)' 
~6~ 77' (w) ^ 8 



du}. 



6(z) 



2 r(f) 
- 2 r(^) 



l + £- 



-1 £2 j^)(3^2_2;2_17) 

~2~ ^ 24 z + 2 

-A+(A-37)^ + 37(^^2)(z + 4) 



From this, residue theorem gives 

— E 



_D-2 



£^0 



r(j + i)Cfl(j + i)S(j) + 



^_^B£^T{D)(:j,{D + 1)A{D) 



27rai (2£) 

1 Wd-D-2 

2-Kax (2£) 



27rai (2£) 

^_2rp-i)CH(i?-i)6(i?-2). 



— ^^rp - 2)CfiP - i)^p - 2) 



duj 



V{z)2-%n{z + 1)£-^ (Ch (z - .? - 1; f ) ^(z) - zS(z)Ch (z - j + 1; f )) 



1 5 

For •perfectly conducting boundary conditions on both spheres, we have ^ ~' ^ ^24 ^'^^ modes and 
D — 2 3 7 

X = — o 1 7 ^ 777 fo'^ TM modes. Using ([T8| . we find that if Z? > 4, the first three leading terms of tlie TE and 

2 . ? 24 

TM contributions to the zero temperature Casimir free energy are given respectively by 



P'-l -1){SD'^ - 2D -17) 2^'*-4i:»3^20i:'2 + 76P»-21Cfi(P'-l) 



24(1? + 2) 



6D{D-l){D + 2) CniD + l) 



Etm,o ^E^f^Q^ 1 + £ 



D-l ^{D-1){3D'^ ~2D-17) ^ D"^ - AD^ - 16D^ + AD + 87 Cr{D - I) 



2A{D + 2) 



6^(1? + 2) (7^-1) Ci?(^ + 1)(' 
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where EJ^^^q and -Exm^o respectively the proximity force approximations for the TE and TM contributions which 
are {D — 2)/{D— 1) and l/(-D — 1) times the total proximity force approximation ([TUl) . When D = 3, the term ■tn/j £i_4 
has to be set to zero. One obtains 

Summing the TE and TM contributions, we find that if I? > 4. the asymptotic expansion of the zero temperature 
Casimir free energy is given by 

p pPFA /.^ , {D-1){3D^-2D~17) , D* - 4ZJ3 + 20i^^ + 40D - 129 Cfl(^ - 1) 1 

£'Cas,o'-i^Ca..oji+e 2 24^ + 2) ^ 6i^p + 2)p-l) Ci?(-D + l)J' 

and if D = 3, 

^--°^-18(k^(^ + ^^+15+^4^ + - 

4- D 3 7 

For infinitely permeable boundary conditions on both spheres, A = — —^-f— — for TE modes and A = 

1 5 

-, 7 = — — for TM modes. Therefore, wc find that ii D > 4, the first three leading terms of the TE and TM 
contributions to the zero temperature Casimir free energy are given respectively by 

Etk,o -E^E.,<l + .^ + e 24(i^ + 2) ' 6D{D + 2) Cr{D + I) [' 



^'TM.o '^-E'TM^ 1 + £ — :z ^ . — £ 



^TM,0 "^-C/TM,0 



D-1 ^{D -1){3D^ -2D -17) ^[D^ - 7){D - 3) Cr{D - I) 



24(L' + 2) QD{D + 2) Cr{D + 1) 



When D = 3, the TE contribution is given by ([27]), and the TM contribution is given by (^5)1 due to duality. Summing 
the TE and TM contributions, we find that if Z3 > 4, the asymptotic expansion of the zero temperature Casimir free 
energy is given by 

pPFA D-l ^{D - 1){3D^ - 2D - 17) ^D"^ - AD^ + 32D^ - %D - %3 C,r{D - l)\^ 

Cas,o~ cas.oj + + ^ 24(i? + 2) ^ QD{D + 2){D-l) Qr{D + 1)]' 

It is interesting to note that in the case of homogeneous boundary conditions, the first analytic correction to the 
Casimir free energy has the form 



E = ijPFA ( 1 + e 



D- 1 



both in the high temperature region and the low temperature region. This is already true for the TE and TM 
contributions separately. It also follows that the proximity force approximation always underestimates the strength 
of the force. 

B. Mixed boundary conditions 

In this case, ai 7^ ol2tI3i 7^ 132- The expression for Au{v{lj)) is more complicated: 

Q2k-i{t{[l + e]u)) - P2k-l{t{uj))\ 



A^{vuj) - exp -2v {ri{[l + £\uj) - r/(w)) - 2 ^ 



fc=i 



The polynomials Pk{t) and Qk{t) are equal to Di^{t) or Mk a{t) depending on the type of modes and the boundary 
conditions, as shown in Table Hill 
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TABLE III: The polynomial Pk{t) and Qk under different boundary conditions. 



Type of mode 


Sphere 1 


Sphere 2 




Qk{t) 


TE 


perfectly conducting 


infinitely permeable 


Dk{t) 


M. i-D(t) 


TM 


perfectly conducting 


infinitely permeable 




Dk{t) 


TE 


infinitely permeable 


perfectly conducting 






TM 


infinitely permeable 


perfectly conducting 


Dk{t) 





In the present case, the computation is more involved because Pk{t) ^ Qk{t)- Proceed as in the previous section, 
we find that 



D-2 



En 



2'Ka 



— E ^^^^-^ r(z)2-^(l - 2-^)Ci?(z + 1)£-M - J - 1; A{z) zC{z)Cn {z-j-r 1; f ) 

1 z(z + 1) 



^2 2 -^WC//(^-J+3;f) 



(28) 



where ^(z) is the same as before, 

/•OO 



77'(u;) 

/•OO 

Q{z) = / (c^77'(c^))-^-2^(t(w))2du;. 

"'0 



2 Vj'iiS)^ £ ujrj'{uj) 



dbJ, 



Here T(t) = Qi{t) — Pi{t) can be written as T{t) — 6t + kI^ . On the other hand, write Qi{t) as Qi{t) = Xt + jt^, we 
find that 



C(z) 



Then 

£^0 



- 2 r(^) 

/ , ^ + 3 2 + 3)(^ + 5) 
' 2 r(^+^) +2dA,— : + K 



z + 4 (z + 4)(z + 6) 



z + 1 \ (z + l)(z + 3) 1 



(5 + K 



2 + 2 



n-2 



^ E ^iS^ro- + 2)(i - + 3)AU + 2) - J- 5: r(, + 1)(1 - + DCO') 



27rai ^ 2J+2eJ+2 

_D-2 



3=3 



27rai 2^£'D 



r(i?)(l-2-^)Cfl(Z? + lM(i?) 



L>;£)-4 



rp - 2)(1 - 22-^)Cflp - 1M(-D - 2) 



27rai 2^-2e^-2 



1 ^D:D-2 



27rai 2^-2e-D-2 



^r(i? - 1)(1 - 2^-^)Cfl,(i? - l)C(i? - 2) 



1 _ 2)(1 - 2^-'')U{D - 3)giD - 4) 



27rai 2^-3e^-2 

If the sphere r = ai is perfectly conducting, and the sphere r — a2 is infinitely permeable, we find that the first three 
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leading terms of the TE and TM contributions to the Casimir free energy are given respectively by 

2 {D - 1){3D^ - 21) - 17) 2 5-0^ - ISD^ - 59D - 15 2^ - 4 (niD - 1) 
24(1? + 2) 6l)(IJT2) 2P~10i(dTT) 

2 2(i:''^ - eL'^ ^ 21)2 281) - 13) 2^ - 16 Ci?,(£' - 3) 1 

+ ^ i?(i?-i)(i?-2)(i) + 2) 2^-1 Ci?(^ + i) r 



^™^° -^^?M^o jl + + e 2^-1ChP + 1) 

2 (D - l)(3i:i2 - 27:> - 17) 2 S-D-'' - SI)^ - 23L» + 9 2^ - 4 CflXl? - 1) 
24(ZJT2) 6i:'(L> + 2) 2^ - 1 Cr{D + 1) 



D-1 2{D^ - 2D - I) 2" - A (^r{D ~ I) 



2{D + 1){D^ - 31)2 _ 3j) _^ 11) 2D - 16 - 3) 
{D-1){D -2)D{D + 2) 2D - 1 Cj^{D + 1) 



where -E^e'^ ^'^'^ ^tm\) ^'"^ respectively the proximity force approximations to the TE and TM contributions which 
are (D — 2)/(D — 1) and 1/{D — 1) times the total proximity force approximation (fTTj) . It follows that the total zero 
temperature Casimir free energy is 



i?-l 2(j)-3) 2^-4Cflp-l) 
2 i) 2^-lCi?P + l) 



2 (D - 1)(3Z)2 - 21) - 17) 2 ^D'^ - SOi?^ - 321)2 _^ ggZ) + 39 2^ - 4 Cfl^i) - 1) 



£ ^ . h £ 



24(Z) + 2) 6D{D-l){D + 2) 2^-lCfl(£' + l) 

2 2(1)'' - 61)3 ^ g£,2 _^ 24Z) - 37) 2^ - 16 Cr{D - 3) 1 
' £'(i:'- 1)(D- 2)(L> + 2) 2^-1 Cb.{D + 1) \ ' 



When D = 3, we have 



77r3 I / 40 \ £2 13 192 

^2880^1^ + - T^j + 15 - ^"t^ +^^V|' 

77r3 / 40 \ £2 27 2 192 1 

^™^° -2880^1 ^ + + 7;? j + 15 + ^ V + ^ V| = (2^) 

77r3 2/^1 1 192 ' 

-^1440^1^ + ^ + ^ U + ^+7^ 

If the sphere r = ai is infinitely permeable, and the sphere r = 02 is perfectly conducting, we find that the first three 
leading terms of the TE and TM contributions to the Casimir free energy are given respectively by 

D-1 2(D2-4Z) + l)2^-4Cfl(D-l) 



^^TE,o "^^TE^^ 1 + ^ — £ 



D{D - 1) 2^-1 Cr{D + 1) 



2 (D - 1)(3Z)2 -2D - 17) 2 7-0"^ - 28Z)3 + 8i:»2 _^ 148^) - 63 2^ - 4 Cfl(D - 1) 
''"^ 24(Z) + 2) ^ 6Z)(Z) - l)(D + 2) 2^ - 1 Ci?(D + 1) 

2 2(1)4 _ Qjj3 ^ 21)2 ^ 28Z) - 13) 2^ - 16 Cr{D - 3) 1 



+ £ 



Z)(L> - 1){D -2){D + 2) 2D -I (j^{D + 1) 
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D-l 2{D^ - 2D - 1) 2^ - A Cr{D - 1) 



■ e 



2 D{D-l) 2D - I Cj^{D + 1) 

[D - l){iD^ -2D - 17) 2 '^D^ - ISi?^ - 40i:'2 + 64D + 57 2^ - 4 (^r{D - 1) 
24(L> + 2) ^ 6l?(I)~"T)(iJT2) 2^ - 1 Qr{D + 1) 

2 2(D + 1){D^ - 3D^ -:iD + 11) 2^ - 16 C^rI^* - 3) 



[D - l)iD - 2)D{D + 2) 2D -I ^j^iD 
It follows that for the total zero temperature Casimir free energy, 




-Ecas,0 ^E?;l^ri{ 1 + e £ 



D-l 2{D - 3) 2^ - A Cr{D - I) 




2D - 1 Cr{D + 1) 

2 ID"^ - 28D^ + 20D'^ + 112D - 183 2^-4 (i^{D - 1) 



■ £' 



24(1) + 2) 6L'(7:>-l)(i:» + 2) 2-° - 1 Cij(L> + 1) 



2(1)4 _ g^3 ^ g^2 24L» - 37) 2D - 16 C,r{D - 3) 



i?(Z? - \){D - 2){D + 2) 2D -I (jj{D + 1) 



The D = 3 ease can be obtamed by duality from ([29]) as before. 

As in the high temperature region, the corrections to the proximity force approximations in the case of mixed 
boundary conditions are more complicated than the case of homogeneous boundary conditions. We find that the 
first correction terms are already different for different combinations of boundary conditions. Again, between the 
two scenarios of mixed boundary conditions, we find that the force is stronger when the smaller sphere is perfectly 
conducting and the larger one is infinitely permeable. 



VI. THE LOW TEMPERATURE ASYMPTOTIC EXPANSION OF THE THERMAL CORRECTION 

To find the low temperature asymptotics of the thermal correction, we use the Abel-Plana summation formula 
[70I - I72I I. which states that for a well-behaved function g{z), 

2^(0) + Y^gip) = g{x)dx + 1 e2^y„i 

+ 2^^ ^(^)Res. ( _2ff_ 1 - 27r» ^i^)^^^A^Mrr\ 

Ro2>0, Im2>0 ^ ^ Ro2>0, lmz<0 ^ 

where w{z) = 1/2 if Re z = and w{z) = 1 if Re z > 0. Applying this to the Casimir free energy ([S|) gives 

E = Eii + — di{D) [ ( ^) ^ — + exponentially decaying terms, 
27r ^ Jo eT - 1 



(30) 



where /; (^) is given by ([4]) . This formula can also be obtained by deforming the contour of integration in ([6]) from the 
positive real axis to the imaginary axis. The middle term in pop is the term that would give the leading terms to the 
thermal correction in the low temperature region. These can be obtained by expanding /i(j^) — fi{—iC) in ascending 
powers of ^ and applying the formula 

di^T{^i + l)CR{^x + l)T^'+\ 



eT — 1 



Observe that the even powers of ^ in f{£) would vanish in the expression fi(iS,) — fi{—i£,)- Therefore, we only need 
to concentrate on the terms in /(^) that is not even in As in l5ll, we use the following smaU z-expansion for Iu{z) 
and K^{z): 

/.W = (|)^r(TT^(i + ^.W), (31) 
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'z\^i+D--2 r(i - v) 



1 / z\~^ 



2 V2 



rH i + ^,.(z) + 2(-i) 



2/+D-2 



1 



l + ^^(z)] , £> odd 



D even. 



(32) 



Here ^u{z) and arc functions that only contain positive even powers of z. From (PT|) and (|32p. it follows that 



where 



1 



0.1 — 0.9 



1 1 ai+Piu 



"2 - fi2V 



ai ~ Pif Q2 + /32'^ / ^2 



2v 



From these, we can see that the leading order term of /z(i^) — is of order ^2;+-D 2_ Therefore, when dT <C 

^ 1, the leading thermal correction comes from the term with / = 1. This implies that 



AtS 



^ ' ' -afr(i:> + i)Cfl(L> + i)r^+i 



CDn2 2^-1 2ai - /Si!) 



DF 



di(jJ)F(^) 2ai+/3ii^ 
V^ai F(f) 2ai-/3ii? 



CHp + i)(ar)^+' + 



Notice that the leading term in the thermal correction does not depend on the boundary conditions and the radius 
of the larger sphere. Using the fact that 



foi(-D) = A 



we find that if the smaller sphere is perfectly conducting, then the leading terms of the thermal corrections of the TE 
contribution, the TM contribution and the total Casimir free energy are given respectively by 

A^i^TE - - 4^^^i^C«(I) + l)(aiT)^+i + . . . , 



T-C/TM 



At£^, 



2Y^ai 
^o, F(f) 
20Fai F(f) 



<fl(i^ + l)(aiT)^+^ + 



(33) 



<i?(i^ + l)(aiT)^+i + 



Note that the leading term of the TM contribution is always negative twice the leading term of the TE contribution. 
If the smaller sphere is infinitely permeable, then 



At^^te 



'(I?-2)V^ai F(f) 



<i?(^ + l)(aiT)^+i + 



AtStm ^^^^r^Ci?P + l)(aiT)^+i + 



/Trai F 
D 



(34) 



(Z?-2)V^ai F(f) 



F f-^^ 

^^Ci^(^ + i)(«ii^)''+' 



When I? > 4, we notice that the thermal correction is larger when the smaller sphere is perfectly conducting. When 
= 3, we have specifically 



3rTi4 



Ar-Ecas ^ —ra\T 
15 



(35) 
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for any boundary conditions. Notice that this is the negative of the low temperature leading term of the thermal 
correction to the Casimir free energy of a single perfectly conducting sphere of radius Oi (Ullll. In fact, one can show 
that (|33p and ((34|) give respectively the negative of the leading thermal correction for a single perfectly conducting 
sphere and a single infinitely permeable sphere in i'-dimensional space. 

In the case of two infinite parallel plates, the low temperature leading term of the thermal correction to the force 
density is given by [67l | : 



-D + l 
TT 2 



-Cr{D + \)T 



D + l 



regardless of the boundary conditions on the plates. Multiplying by the area of the sphere of radius oi, the proximity 
force approximation of the low temperature leading term of the thermal correction is 



' V^al r(f) 



<fl(i^ + l)(air) 



D+l 



(36) 



However, the leading terms (|33p and ([M)) derived from the exact formulas of the Casimir free energies showed that 
when the smaller sphere is perfectly conducting, then the low temperature leading term of the thermal correction to 
the force is 



r(#) 



■Ci?(i? + l)(aiT) 



D+l 



(37) 



Whereas if the smaller sphere is infinitely permeable, then the low temperature leading term of the thermal correction 
to the force is 



r>2 rf-^) 



{D-2)^al r( 



(38) 



As expected, these do not agree with the temperature correction in the proximity force approximation p6p . The point 
is that under the low temperature condition dT <C <^ 1, the thermal correction is much smaller than the zero 
temperature Casimir energy by an order and, thus, the proximity force approximation is not applicable. However, 
the proximity force approximation remains applicable for the calculation of the total free energy as we have seen in 
Section |Vl Similar situation has been observed in D = 3 dimensions for the case of a sphere in front of a plate 51 j . 



VII. 



CONCLUSION 



In this article, we studied the Casimir interaction between two concentric spheres in {D + l)-dimensional spacetime 
due to the confinement of the electromagnetic field between the spheres. We consider the cases of perfectly conducting 
- perfectly conducting, infinitely permeable - infinitely permeable, perfectly conducting - infinitely permeable and 
infinitely permeable - perfectly conducting boundary conditions on the spheres. The first two are referred to as 
homogeneous boundary conditions, and the last two are called mixed boundary conditions. For homogeneous boundary 
conditions, the Casimir interaction between the spheres is always attractive. For mixed boundary conditions, it is 
always repulsive. 

We arc particularly interested in studying the asymptotic behaviors of the Casimir free energy when e, the ratio of 
the separation between the spheres to the radius of the smaller sphere, is small. Both the high temperature region and 
the low temperature region are considered. In the high temperature region, the Casimir free energy is dominated by 
the classical term which is the term corresponding to the zeroth Matsubara frequency. In the case of two concentric 
spheres, this term is quite simple. It can be written as a series in elementary functions. In the low temperature 
region, the Casimir free energy is dominated by the zero temperature term, which has to be expressed in terms of 
Bessel functions. The first three leading terms are computed explicitly. In the high temperature region, the leading 
terms are of order Te^~^, and they coincide with that obtained using proximity force approximation. For the zero 
temperature terms, the leading terms are of order e~^, and they also agree with the proximity force approximations. 
It is interesting to observe that the asymptotic expansions of the Casimir free energies have the following universal 
structure: 



Cas 



PFA 
Cas 



1 



D-1 



-£^o(-D) 



C{D + i-2) 



'^2(-D)+"^i(-D) 



C{D + i-2) 



C{D + i- 4) 



C(^ + *) 
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where i = I'm the low temperature region, and i = in the high temperature region. ^Wq{D),^q{D), '^i(D), '^2{D) are 
rational functions of D that are 0{D'~'), 0{D^), 0{D), 0{D^) when D is large. In the case of homogeneous boundary 
conditions, the terms 3§o{D) and 'r^o{D) are absent. 

In general, the corrections to the proximity force approximations are more complicated in the case of mixed boundary 
conditions compared to the case of homogeneous boundary conditions. In fact, for homogeneous boundary conditions, 
the first correction is the same when the two spheres are both perfectly conducting or both infinitely permeable. For 
the two scenarios of mixed boundary conditions, the first corrections are different. It is observed that the Casimir 
interaction is stronger when the smaller sphere is perfectly conducting and the larger sphere is infinitely permeable. 

Finally, the low temperature leading terms of the thermal corrections to the Casimir free energies are computed. 
They are finite when e — >■ 0+ and are of order T^+^. It is interesting to find that these leading terms are independent 
of the larger sphere. They do not depend on the radius or the boundary conditions on the larger sphere. As have been 
observed by a few researchers for other geometric configurations [5l|, [Sg , this case is outside the application region of 
the proximity force approximation. 

For future works, it would be interesting to consider eccentric spheres and compare the results with those obtained 
here. 
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